Metamagnetic phase transition in the Ising plus Dzyaloshinskii-Moriya model 
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We study the ID ferromagnetic Ising (spin- 1/2) model with the Dzyaloshinskii-Moriya (DM) 
interaction. We analyze the low energy excitation spectrum and the ground state magnetic phase 
diagram using the Lanczos method. The DM interaction-dependency is calculated for the low-energy 
excitation spectrum, spiral order parameter and spin-spin correlation functions. We show that a 
metamagnetic quantum phase transition occurs between the ferromagnetic and spiral phases. The 
existence of the metamagnetic phase transition is confirmed, using the variational matrix product 
states approach. 

PACS numbers: 75.10.Jm Quantized spin modcls;75.10.Pq Spin chain models 
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I. INTRODUCTION 



Studying one dimensional quantum spin systems have 
been obtained many interesting results. The Ising spin 
models pose intriguing theoretical problems because an- 
tifcrromagnctic (AF) and ferromagnetic (FM) systems 
with spin-1/2 have a gap in the excitation spectrum. 
Therefore they reveal an extremely rich behavior dom- 
inated by quantum effects. In particular, the spin-1/2 
Ising model in a transverse magnetic field (ITF) dis- 
plays a pragmatic example of a quantum phase transi- 
tion. Theoretically, the ITF problem is exactly solved^ 
and found a phase transition at a finite value of the trans- 
verse magnetic field: h c . This is a quantum critical point 
and the phase transition is continuous in nature. 

The antisymmetric spin exchange interactions between 
spins, known as the Dzyaloshinskii-Moriya (DM) interac- 
tion, play an important role in physics of spin systems 2 ^. 
The DM interaction idea originated from the deviation of 
experimental data from the theoretical predictions, based 
on the Heisenbcrg spin Hamiltoniana^ ' 6 ' 7 ' 8 ' 9 ' 10 ' 11 ' 12 . 
Generally the DM interaction between two spins Si , S2 
can be written as D- (Si x S2) with an axial DM vector D. 
In actual systems, the direction of D vector is fixed by 
the microscopic arrangement of atoms and orbitals. In a 
spin chain, the DM vector may vary both in direction and 
magnitude. However, the symmetry arguments usually 
rule out most of possibilities and confine the theoretical 
discussion to two principal cases. The first one is the uni- 
form DM interaction, D = constant over the system^ 3 - 
and the second case is the staggered DM interaction^ 4 - 
with anti-parallel D on adjacent bonds. Since the DM 
interaction is rather difficult to handle analytically, the 
effect of this interaction are only partially studied so far. 
In this sense we study the Ising chain (FM) with DM in- 
teraction, that its Hamiltonian (by considering a periodic 



chain of N sites) is given by 



H 



(S j x S j+1 ) 



(1) 



where Sj is spin-1/2 operator on the j-th site, and J > 
(J < 0) denotes the AF (FM) coupling constant. In 
a very recent work, using the quantum renormalization 
group and numerical Lanczos methods, the ground state 
phase diagram of the AF Ising chain ( J > 0) is studied^. 
It is shown that the ground state phase diagram consists 
of AF and spiral phases. 

By considering uniform DM vector as D = Dz, and 
doing the rotation about z axis as S^ — ► Sj t e ±i 5 !L , the 
Hamiltonian is transformed to XXZ chai n 16 ! 17 with an 
anisotropy parameter jy, i.e., 



H 



( QX QX 



qy qy i 



^ QZ QZ 



(2) 



The XXZ chain model was solved by Bethel, and its 
ground state phase diagram is well knowni^. The Neel 
regime is governed by > 1 and there is a gap in the 
excitation spectrum. For jy < — 1, the ground state is in 
the FM phase and there is a gap over the FM state. In 
the region — 1 < jy < 1, the ground state of the system 
is in the gapless Luttinger liquid (LL) phase. Thus, by 
increasing the D value, the system undergoes a quantum 
phase transition from the gapped Neel (or FM) phase to 
the gapless LL phase at the critical value D c = J. 

In this paper we study an Ising spin-1/2 chain with 
FM exchange (J < 0) and the uniform DM interaction 
using the numerical and analytical approaches. In the 
forthcoming section, we apply the Lanczos method to 
diagonalize numerically finite chain systems with lengths 
N = 8, 10, 24. Using the exact diagonalization results, 
we calculate the spin gap, magnetization and spin-spin 
correlations as a function of the DM interaction. We 
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FIG. 1: Difference between the energy of the two lowest levels 
and the ground state energy as a function of the DM vector, 
for chains with FM exchange (J = —1.0) including different 
lengths N = 12, 16, 20. 

have also calculated the spirality in the ground state of 
the system. Based on the exact diagonalization results we 
obtain the ground-state magnetic phase diagram of the 
model. By taking the DM interaction as the control pa- 
rameter we show that a metamagnctic phase transition, 
can be observed in the ID ferromagnetic Ising model. In 
section III, the observed metamagnetic phase transition 
from the numerical lanczos method is confirmed by the 
analyzing the results of the variational matrix product 
states approach. Finally, we conclude and summarize 
our results in section IV. 



II. NUMERICAL RESULTS 

In this section, to explore the nature of the spectrum 
and the quantum phase transition, we use the Lanczos 
method to diagonalizc chains with length up to N = 24. 
We have computed the three lowest energy eigenvalues of 
chains with FM exchange J = —1.0 and different values 
of the DM vector. To get the energies of the few lowest 
eigenstates we consider chains with periodic boundary 
conditions. In Fig[TJ we have presented results of these 
calculations for the chain sizes N = 12, 16, 20. We define 
the excitation gap as a gap to the first excited state. It 
exhibits that the energy spectrum is gapped at D = 0, 
while by increasing the D the energy gap decreases lin- 
early and vanishes at D c = \ J\ = 1.0. We got this critical 
value as an exact value since there was no finite size cor- 
rection. In the region D < D c , the difference between 
the energy of the first excited state and the ground state 
energy shows an universal linear decrease with increas- 
ing DM vector, which is independent on the chain length 
(within the used numerical accuracy). By increasing the 
DM vector for D > D c , the gap opens again in finite 
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FIG. 2: a) The spin structure factor S zz (q) as a function of q 
for FM Ising chain with DM interaction (J = —1.0 and N = 
24). b) The spin structure factor S zz (q = 0) as a function of 
the DM interaction D for FM Ising chain (J = —1.0) including 
different chain lengths N = 16, 20, 24. 

chains, but by extrapolating finite size results to N — > oo, 
we found that there is no gap in the spectrum. That is 
in agreement with the transformed Hamiltonian results 
(Eq. O. Hence there are two gapped and gapless phases 
in the ground state phase diagram. 

To study the magnetic order of the ground state of 
the system, we have implemented the Lanczos algorithm 
on finite chains to calculate the lowest eigenstate. The 
symmetry breaking cannot occur in finite size systems, 
thus instead of the magnetization, M z = ^2j(Sj), we 
focus on the spin-spin correlation functions. The static 
spin structure factor at momentum q is defined as 

SZZ (l) = ^T, eZqn ^ S Un)- (3) 

n 

It is known that the spin structure factors give us deeper 
insight into the characteristics of the ground state. The 
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FIG. 3: The numerical results on the ground state energy 
(normalized by N) of a FM Ising model with DM interaction 
as a function of D for different chain lengths N = 12, 16, 20 
(|J| = 1.0). 



-1 



i i i i i i i i i i_ 



0.05 0.075 



l/N 



N=16 
N=20 
N=24 

_L 



1.5 
D 



FIG. 4: The spiral order parameter as a function of the DM 
vector for the FM Ising chain (J — —1.0). Inset shows the 
spiral correlation function, C z , as a function of the l/N for 
different values of the DM vector: D = 1.1, 1.2, 1.3, 1.4, 1.5. 



D-dependency of the spin structure factor, S zz (q = 0), 
is qualitatively the same as the uniform magnetization, 
M z . 

In Figf2Ia) we have plotted the spin structure factor, 
S zz (q), as a function of momentum q for a chain with 
length N = 24 and different values of the DM vector 
in the region D < D c . We found that in this region 
the function S zz (q) is independent of the value of D. It 
shows that S zz {q) has a very sharp peak at q = 0, cor- 
responds to the FM ordering. On the other hand, there 
is not any peak in the region D > D c which confirms 
that all spins should aline in the xy plane. In Figf^b), 



we have plotted S zz (q = 0) as a function of D for the 
chain lengths N = 16, 20, 24. For D < D c the spin struc- 
ture factor, S zz (q = 0), is equal to the value 0.25, which 
shows that the ground state of the system is in the fully 
polarized FM phase. This quantity is independent of the 
chains size, therefore the value of the magnetization does 
not change in the thermodynamic limit N — ► oo, and 
the FM ordering is a true long range order in the re- 
gion D < D c . One of the most interesting predictions 
of this model is that the magnetization as a function of 
the DM vector displays a jump for certain parameters. It 
means that the spontaneous magnetization, M z , remains 
at the saturation value in the region D < D c . But at the 
critical value D = D c , the spontaneous magnetization 
jumps to zero. A rapid increase (or discontinuity) at crit- 
ical value of the control parameter in the magnetization 
curve is called the metamagnetic phase transitio n 20 ' 21 . 
This phenomena that observed in the ground state phase 
diagram of the ID frustrated FM Heisenberg mode l 22 ' 23 , 
has defined the phase transition between the FM and AF 
phases. Surprisingly, however, we find that the metam- 
agnetic phase transition can be observed between the FM 
and spiral phases in the ID ferromagnetic spin-1/2 Ising 
model with DM interaction. In following we draw a sim- 
ple physical picture for this phenomena. 

The wave function of the ground state in the absence 
of the DM interaction has a form \ipcs) = I TT ■••)• Ap- 
plying a uniform DM interaction on the FM state \ipas) 
yields 



^D.(S j xS j+1 )|te)=0, 



(4) 



which shows that in the presence of a DM interaction the 
ground state is fully polarized and does not change up to 
the critical value, D c . The system is fully FM in < D < 
D c , and lies in the subspace Stot — N/2 with two-times 
degeneracy. But it becomes an incommensurate state for 
D > D c , where the energy gap is strongly suppressed. 
At the critical point, two distinct configurations with the 
energy Eqs = — \N\J\ are the ground states, where one 
is fully polarized in z direction. In Fig[3J we have plotted 
the Lanczos results on the ground state energy per sites 
as a function of D for the chain lengths TV = 12, 16, 20. 
As we expected in the region of D < D c the energy per 
spin of the ground state has the constant value —0.25 and 
independent of the DM interaction (Eqs/N = — |J|/4). 

On the other hand, we have also computed numerically 
the transverse spin structure factors (S xx (q),S vv (q)) for 
different values of DM vector. Due to symmetry S xx (q) 
is the same as S vy (q). We found that for the values 
of the DM vector D < D c , S xx (q) = S yv (q) = in 
well agreement with the saturated ferromagnetic phase 
in the z direction. In the region D > D c , the transverse 
spin structures showed two peaks at q — f,^-, which 
is a justification of the spiral order—. It is showed that 
the DM interaction can induce the spiral phase in the 
ground state phase diagram of the spin systems 2 ^, which 
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is characterized by the nonzero value of the spirality 

x* = ^E<( s * xS '-+i)*>- (5) 

i 

One should note that there are two different types of the 
spiral ordered phases, gapped and gaples o 25 ' 26 . There- 
fore, the definition of the spiral correlation function as 



TV 
n=l 



(6) 



provides further insight into the nature of different 
phases. In FigfJ] we have presented results of these 
calculations for Ising chains with different lengths N = 
16, 20, 24. In complete agreement on the results of mag- 
netization, the spirality in the FM Ising chain shows 
a plateau in zero value at D < D c . Which confirms 
that there is no spiral long range order in the mentioned 
regime and the ground state is in the saturated FM phase. 
The spirality remains close to zero value in the region 
D < D c , but at the critical point, it jumps to a non-zero 
value. Which is also an indication of the metamagnetic 
quantum phase transition that occurs only in the case of 
FM Ising chains. The inset of the FigfJ] shows the spi- 
ral correlation function, C z , as a function of the 1/N for 
different values of the DM vector. It shows that the spi- 
rality remains non zero in the thermodynamic limit for 
D > D c , that corresponds to the spiral long range order 
in the xy plane. 



III. VARIATIONAL MATRIX PRODUCT 
STATES APPROACH 



The matrix product state is defined a: 
*) =Tr(gx ■g2---9N), 



,27,28 



(7) 



where the elementary matrix gj represents the matrix- 
state of the jth spin cell. The size of these elemen- 
tary matrices depends on the problem. In this pa- 
per the simple one dimensional case is considered, i.e., 



fJj 



t)i + bj\ Where a,j and bj are probability 



amplitude for two configuration of the spin at site j, and 
N is the lattice site number. The goal is to determin the 
ground state energy of the FM Ising spin system with 
DM interaction. In this respect, the variational energy is 
obtained by 



E var = (H) = 



{V\H\V) 



E 



H 



[Gj ■ G j+1 ] 



(8) 



where (*|*) = g\ 
defined by Gj = 



® 5i ' ' ' 9n ® 9n 

„ |2 i 



UjGj, and Gj 



is 



® gj 



D • Sj x S k , and Sf = g) 



S?9j- 



Here H 



The minimum of the variational energy function cor- 
responds to the ground state energy of the system. Us- 
ing the normalization condition ( V I'| V I') = 1, variational 




FIG. 5: The ground state energy (normalized by N) of FM 
Ising with DM interaction using the variational matrix prod- 
uct state. Inset shows the spontaneous magnetization \M Z \ as 
a function of the DM interaction, using the variational matrix 
product state (|J| = 1.0). 



parameters can be mapped to a,j = cos(8j)e lifij , bj = 
sm(9j)e lip i . Therefore one can obtain S| = cos(20j)/2, 

S+ = (S-)* = l/2sin(20 J -)e^, where fa = - tpj). 
By choosing D = Dz, it is easily found that 

E var = - Y^[Jcos(28j) cos(2%i) + 

j 

D sm(26j) sm(26 j+ i) sm((/)j - <f> j+1 )]. (9) 

By minimizing the above equation, the ground state en- 
ergy (Egs) m the case J < shall be achieved. One 
can show that the ground state energy has the constant 
value, i.e. Egs = — iV|J|/4 for D < \J\, and it behaves 
almost linearly for D > \J\ as Egs = —ND/A (please 
see the Fig [5]). 

The magnetization can be written as M z = 
Y^jCOs(29j)/(2N), thus by considering the minimized 
variational parameters it follows the metamagnetic phase 
transition curve which has shown in the inset of the FigjSl 
where \M Z \ = 0.5 at D < |J| and zero for the D > |J|. 
Also the spirality is given by 



X = 



± E M%>j) sm(20 j+1 ) sin(0, - <f> j+1 ), (10) 



where using above results, one can obtain \ z 
D < \J\, and \ X Z \ = 0.25 for D > \ J\. 



for the 



IV. SUMMARY AND CONCLUSION 

In this paper the elementary excitations and the mag- 
netic ground state phase diagram of the ID spin- 1/2 FM 
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Ising model with the Dzyaloshinskii-Moriya interaction 
have been thoroughly investigated by numerical tools and 
variational schemes. Using the analytical and numeri- 
cal approaches, we have shown that there are two dif- 
ferent phases in the zero-temperature phase diagram of 
the model. To provide a physical picture of the ground 
state phase diagram of the model, by a redefinition of 
the spin variables the model is mapped onto an XXZ 
Heisenbcrg chain. Where the anisotropy parameter is re- 
lated to the DM vector, and identified a commensurate- 
incommensurate (C-IC) quantum phase transition be- 
tween the gapped and gapless phases. The numerical 
experiment with high accuracy, has shown that in the 
ground state phase diagram of the FM chain with DM 
interaction, there is only one fully polarized FM phase 
below the critical value: D c = \J\. However at the criti- 
cal value, a metamagnetic phase transition occurs to the 
spiral gapless phase. Moreover, we used a variational 
method for evaluating ground state energy. By consider- 
ing ID simple matrix product state, the zero temperature 
energy diagram of the ferromagnetic Ising chain with DM 
interaction has been obtained. Our data shows the sud- 



den jumping in the magnetization and spirality curves at 
D c = \J\. 

Based on the results we conclude that properties of the 
FM and AF Ising spin- 1/2 chains with the uniform DM 
interaction are very different. In principle, the ground 
state phase diagram of the FM chain consists of two 
phases which are gapped and gapless. In the case of 
FM chains, for D < D c , there is only one saturated FM 
phase and at the critical value, D Cl a metamagnetic phase 
transition occurs to the spiral phase. 
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